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Abstract. It has been shown previously that the DM in galactic halos can be explained by a 
self-gravitating system of massive keV fermions (‘inos’) in thermodynamic equilibrium, and 
predicted the existence of a denser quantum core of inos towards the center of galaxies. In 
this article we show that the inclusion of self-interactions among the inos, modeled within a 
relativistic mean-field-theory approach, allows the quantum core to become massive and com¬ 
pact enough to explain the dynamics of the S-cluster stars closest to the Milky Way’s galactic 
center. The application of this model to other galaxies such as large elliptical harboring mas¬ 
sive central dark objects of ~ 10 9 Mq is also investigated. We identify these interacting inos 
with sterile right-handed neutrinos pertaining to minimal extensions of the Standard Model, 
and calculate the corresponding total cross-section er within an electroweak-like formalism to 
be compared with other observationally inferred cross-section estimates. The coincidence of 
an ino mass range of few tens of keV derived here only from the galactic structure, with the 
range obtained independently from other astrophysical and cosmological constraints, points 
towards an important role of the right-handed neutrinos in the cosmic structure. 
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1 Introduction 

The Cold Dark Matter (CDM) model of the Universe, characterized by ordinary matter 
(about 5%), a vacuum dark energy (more than 68%), and Dark Matter (DM, 27%), with an 
equation of state resembling a positive-cosmological-constant (A) type fluid (a.k.a. ACDM 
model) seems to be, at least currently, the cosmological scenario that best fits the plethora of 
the available cosmological and astrophysical data [1]. At present, the nature of DM still elude 
us. Supersymmetry, which provides leading candidates for cold DM, has not been discovered 
as yet, thus prompting us to consider alternative candidates for DM such as axions, or sterile 
right-handed neutrinos with masses higher than 100 keV. 

On the other hand, right-handed neutrinos with masses less than 50 keV may still play 
a role in particle physics today, as conjectured in the so-called right-handed neutrino minimal 
(non-supersymmetric) extension of the standard model (z/MSM) proposed in [2], which has 
been proposed as a viable model for the so-called warm DM (WDM). This model involves 
three right-handed neutrino states, in addition to the three left-handed active neutrinos of 
the standard model (SM) sector, of which the lightest, of mass at most a few tens of keV, 
can live longer than the age of the Universe, thus constituting a viable DM candidate. Such 
relatively light right-handed neutrinos appear compatible with cosmological DM and Big- 
Bang-Nucleosynthesis constraints, provided their mixing angles with the active neutrinos of 
the SM sector are sufficiently small, as shown in figure 1. In general terms, the model appears 
to be consistent with a plethora of diverse astrophysical and cosmological data [3-6]. 

From the astrophysical point of view, the long-pursued study of galactic DM within the 
context of fundamental physical principles including thermodynamics and statistical physics 
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Figure 1. Cosmological constraints on the mass (Mi) and mixing ($i) parameters of the light¬ 
est sterile neutrino state Ni of the uMSM model, consistent with all the current astrophysical and 
cosmological data [2-6]. This picture was taken from the original version in [4]. 

dates from several decades already (see e.g. [7]), given that DM halos show clear universal 
properties ([8] and refs, therein) and are well fitted by different phenomenological profiles 
that resemble isothermal equilibrium spheres [9-12], Due to the collisionless nature of DM 
particles at large scales, it has been recognized that the main mechanisms for the (quasi) 
relaxation of a DM halo within observable time-scales are collisionless processes such as 
phase-mixing and violent relaxation [7, 13]. In contrast to the standard collisional scenarios 
relevant for (stellar-dominant) globular clusters, violent relaxation takes place within a much 
shorter timescale, appropriate for the time-varying (global) gravitational potential, leading 
the system to a long-lived quasi-stationary-state (QSS) that, as shown in [13-17], under 
sufficient mixing conditions can be well described in terms of Fermi-Dirac statistics. More 
specifically, this kind of Fermi-Dirac (coarse-grained) phase-space distribution can be ob¬ 
tained from a maximization entropy principle at fixed halo mass and temperature as clearly 
demonstrated in [14-16, 18, 19], 1 , leading to mass-density equilibrium distributions typically 
composed by a condensed core surrounded by a dilute halo. This mechanism was first de¬ 
rived for classical particles (i.e. stars) in [13, 23], and then extended for indistinguishable 
fermionic particles in [15, 24], More recently, and within the context of fermionic DM, it has 
been argued [11, 12, 25-30] that a system of self-gravitating fermions, which we have referred 
to as inos, with masses in the keV regime, plays an important role in galactic structures. 
In the more general case of fermionic models allowing for central degeneracy [12, 28-31], 
the density of the inos, which we propose here to be identified with right-handed Majorana 
neutrinos 2 , shows three physical regimes as a function of the distance from the center: 1) 

x This procedure certainly implies the necessity for these QSS to be bounded in radius. This kind of 
configurations can be easily obtained within our model, by setting a cutoff in the momentum space of the 
original Fermi distribution as first realized in [20], or, more recently and within our relativistic formalism in 
[21, 22]. In any case the main conclusions presented here do not depend on such a cutoff that only fixes the 
outermost halo boundary conditions. Therefore, we shall adopt throughout this paper (for simplicity) the 
standard Fermi-Dirac statistics, with the momentum cutoff set to infinity. 

2 These neutrinos could be of the DM type appearing in izMSM [2], but such an identification is not binding. 
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an homogeneous inner core where the inos are in a degenerate state; 2) the uniform density 
of the core is followed by a steeply decreasing density and an extended plateau in an inter¬ 
mediate region in which the ino’s description needs still some quantum corrections; 3) and 
finally it ends with an asymptotic p oc r^ 2 classical Boltzmann regime. The latter regime 
is the responsible for the flatness of the rotation curves and therefore it has to fulfill, as an 
eigenvalue problem, a defined value of the circular velocity. It was further shown in [30] that 
this eigenvalue problem allows to determine the mass of the ino as well as the radius and mass 
of the inner quantum core. This kind of core-halo structure for DM in galaxies is consistent 
with the results obtained in [14, 15, 18] within a pure statistical approach. Interestingly, 
similar core-halo distributions have been obtained in modern 3-D numerical simulations in 
the framework of quantum-wave DM approaches [33]. Moreover, such structures appear to 
characterize more generally long-range collisionless interacting systems, including plasmas 
and kinetic spin models [34]. 

The initial conditions for such a core-halo galactic QSS are provided by the aforemen¬ 
tioned collisionless relaxation processes, which imply specific fermionic phase-space distribu¬ 
tions as the ones used here. These quantum fermionic distribution functions, with relatively 
large values of central degeneracy parameters, stabilize the galaxy structures by avoiding a 
thermodynamic runaway, and thus the gravothermal catastrophe, thanks to the Pauli prin¬ 
ciple. This is in contrast to what happens in the case of Boltzmann-like configurations [11], 
where a gravothermal catastrophe (similar to one occurring in globular clusters) is an in¬ 
evitable outcome, even for collisionless DM particles [18, 35]. The precise choice of the free 
parameters for the appropriate Fermi distribution, once the QSS is achieved, is dictated by 
the correct values of the observed halo parameters (circular and/or dispersion velocities, to¬ 
tal mass), as well as the desired quantum core mass and compactness, to account for the 
observed properties of the central massive object. 

The novel approach introduced in [30] was applied to different types of galaxies ranging 
from dwarfs to large spiral galaxies. For ino masses m ~ 10 keV/c 2 , one finds excellent 
agreement with the DM halo observables (see [12, 30], for details). At the same time, the 
approach is capable of providing a theoretical correlation between the inner quantum core 
and the halo mass, which can be compared with observations [30]. We also evaluated the 
possibility of an alternative interpretation to the black hole in SgrA*, in terms of the high 
concentration of DM in the inner quantum core. We concluded that, although a compact 
degenerate core mass M c ~4x 10 6 Mq is definitely possible with an ino of m ~ 10 keV/c 2 , 
the core radius is larger by a factor ~ 10 2 than the one obtained from the observational 
limits imposed by the of S-star trajectories such as SI and S2 orbiting around SgrA* [36, 37]. 
To solve this problem, we propose here the inclusion of specific (self) interactions among 
the inos, which, as we shall demonstrate below, allows for higher central degeneracies and 
higher compactness of the inner quantum core. At this point it is important to stress that, 
already in [23], two-particle interactions were predicted to be non-negligible within the dense 
degenerate cores, due to the appearance of the exclusion principle, in agreement with the 
ansatz considered here. Moreover, the necessity for considering self-interactions in dense and 
very-low temperature fermionic systems, such as the ones studied in the current work, has 
been proven in laboratory experiments. Indeed, in [38], it was argued that the behavior of 
ultra-cold atomic collisions in (effective) Fermi gases, such as 6 Li, can be explained in terms 

Indeed, our inos can be also identified with sterile neutrinos which do not rely on active-sterile mixing, as the 
ones analyzed in [32], and thus consistent with all current cosmological/astrophysical constraints for masses 
in the keV - MeV range, similarly to the range obtained in section (3.2). 
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of a grand-canonical many-body Hamiltonian with a term accounting for the (spin-enhanced) 
fermion-fermion interaction. At temperatures of a fraction of the Fermi energy, or, equiva¬ 
lently, for thermal de-Broglie wavelengths larger than the inter-particle mean distances, the 
self-interactions of the fermions become relevant, in analogy with the situation encountered 
in our self-interacting neutrino model. However, while in the case of laboratory physics an 
external trapping potential (such as the one due to magnetic fields) is needed, in the context 
of DM in galaxies, trapping is ensured by gravity. 

The idea of self-interacting DM was first implemented in [39, 40] for cold DM particles 
with rest masses above 1 MeV/c 2 (and up to 10 GeV/c 2 ), consistent with the nature of 
the interactions and the mean free paths considered. In those works, self-interactions were 
applied uniquely at DM halo scales with typical densities of order 10 _2 Mq/pc 3 , suggesting 
that total cross-sections over the particle mass of order a/m ~ 0.1 — 100 cm 2 /g, would imply 
observational effects in the inner halo regions. It was further shown that a self-interacting 
DM regime with these values of a/m, would generate shallower inner DM profiles, with a 
consequent reduction in the amount of sub-structures, thereby alleviating important prob¬ 
lems of collisionless ACDM simulations, such as the core-cusp [41] and the missing satellite 
problems [42]. However, at the same time, some tension with upper limits in the DM cross 
sections obtained from lensing analysis at galactic cluster scales has emerged. More recently, 
in [43], motivated by a more refined analysis of the Bullet Cluster [44], a set of cosmological 
simulations within CDM were performed, with the aim of studying further the effects of 
self-interacting dark matter (SIDM) on density cores of galaxies and galaxy clusters, con¬ 
cluding that a/m ~ 0.2 barn GeV” 1 = 0.1 cm 2 g -1 is consistent with all the observational 
constraints. 

In the above works, the interactions of DM were modelled by pure classical mechanics 
descriptions, without making any reference to the details of the interactions. By the con¬ 
trary, in the present paper, we analyze the possible consequences caused by a self-interacting 
relativistic field theoretical model of Majorana fermions, with vector type interactions and 
fermion rest-masses in the keV/c 2 range, which can play the role of WDM in galaxies. In par¬ 
ticular, we maintain the collisionless nature of the DM fermions at halo scales, and study the 
two-particle self-interaction effects for different interaction strengths, but only in the (sub- 
pc) region, where the dense fermionic quantum core arises [30], reaching central densities as 
large as 10 16-23 Mq pc” 3 . Our method is fundamentally different with respect to the one in 
other approaches (see, e.g., [43]). In such works the ansatz of self-interacting DM is assumed 
to study its effects on the central parts of halo density profiles coming from cosmological 
simulations, while we use that ansatz to study the self-interacting effects on the quantum 
cores of the DM profiles arising from first principle physics such as quantum statistics, ther¬ 
modynamics and gravitation. At halo-distance scales, our assumption of a collisionless DM 
is understood by the fact that, the non-interacting fermionic DM distribution, leads natu¬ 
rally to cored inner halos (see figures 2-4 below). Thus, within our fermionic DM model 
there is no need to make use of the SIDM hypothesis for the halo, as needed in standard 
ACDM cosmological simulations to alleviate the tension with observations. At the end of the 
concluding section 4, we shall comment on how the above mentioned core-cusp discrepancy 
may be tackled by the current model, and compare briefly our results with predictions made 
by the standard cosmological N-body simulations of cold and warm DM models. We should 
stress though that, although many features related to the galactic structure, where the CDM 
model was challenged, may be explained by our self-interacting fermion (right-hand neutrino) 
model, the latter should not be viewed as a panacea for solving all the current open issues 
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regarding that front. There may well be more than one DM species in the universe, and 
in this work we offer a proposal towards a solution to some important problems in galactic 
structure within the context of our self interacting (right-handed neutr)- “inos”. 

Both approaches, those based on standard cosmological simulations and ours, can be 
thought as complementary in attacking the problem of the distribution of DM in galaxies. 
Thus, it is of interest to compare and contrast any theoretical prediction of a self-interacting 
DM model with the parameters inferred from observations (and/or simulations), such as 
the total cross-section per unit mass a/m . For this purpose, we consider in this work DM 
self-interactions mediated by a massive-vector in the dark sector of minimal extensions of 
the Standard Model (SM), with a SU (3)xSU (2) l.rxU (1) group invariance, and compute the 
total cross-section a through an electroweak-like formalism ( see, appendix B, for details). 
Then, we compare and contrast our theoretical predictions with the observational constraints 
in section 3.4. We show that, the requirement that the cross-section per unit mass agree 
with the one constrained from observations and SIDM cosmological simulations [43], implies 
an allowed interaction-strength window Cy = (gv/ m v ) 2 € (2.6 x 10 s , 7 x 10 8 )Gi? (where 
Gf ~ 10~ 5 GeV -2 is the Fermi coupling constant of the weak interaction), for particle 
masses in the range m € (47,350) keV, which as we argue here are in agreement with 
Galactic-core observables. Here gy and my are the coupling constant of the interaction and 
the mass of the vector-meson mediator, respectively. Since our approach is not bound to 
standard ACDM-based conclusions, we also discuss further constraints on the interaction 
coupling associated with more extreme quantum-core effects. By linking the cross-section 
with the scattering-rate per particle T, we also show here an absolute lower bound for the 
DM interaction strength (and a), by calculating the scattering probability among the inos 
to occur at least once in the age of the galaxy. In this way, we obtain a minimum value for 
the coupling constant Cy. 3 * 

The structure of the article is as follows: in the next section 2 we introduce the model 
of right-handed (Majorana) neutrinos with vector self-interactions, which could be either due 
to a vector field or describe contact current-current type of interactions. A numerical study 
of the induced core-halo structure for galaxies, assuming that the above model is the correct 
one to describe the DM in the Galaxy, is given in section 3. The effects of the self-interaction 
in ensuring higher central degeneracies and higher compactness of the inner quantum galactic 
cores are demonstrated. Finally, discussion of the results and outlook are presented in section 
4. In particular, we specify the order of the vector interactions field strength as well as the 
minimum value of the fermion masses (~ 47 keV/c 2 ) for the model to provide a description 
of the core-halo structure in a variety of galaxies, from spiral to large elliptical, in agreement 
with observations. In case the inos are identified with the lightest right-handed Majorana 
neutrino in the zzMSM, then there is only a narrow (but non trivial) regime of masses for 
which the model can be consistent with astrophysical/cosmological/galactic data in the sense 
considered in this paper and in particle physics applications of the i^MSM. Some technical 
aspects are given in appendices. 


3 It is understood that the coupling CV “runs” with the energy of the SIDM particle, but here we give 

its value technically at zero momentum. Within our low energy approximations the running of CV with the 

energy is very soft and negligible. 
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2 Self-interacting right-handed neutrinos 

We consider in this work a model for SIDM that is a minimal (non supersymmetric) extension 
of the Standard Model with a sterile neutrino. Such models are reminiscent (but different) of 
the general idea of the I'MSM [2-4], Unlike zzMSM, we allow our right-handed neutrinos to 
be self-interacting. In particular we concentrate on the lightest of the right-handed Majorana 
neutrino N\. which plays the role of DM, and we introduce phenomenologically, self-neutrino 
interactions through a massive-vector-meson mediator. Our results regarding the DM 
particle creation mechanism, are only subject to the assumption that the fermions are of 
Majorana type (but the formalism is readily extendable to Dirac). This is the common 
feature we share with the i/MSM, which from our point of view is another interesting and well- 
studied case of sterile neutrinos that we use for comparison, given the intriguing similarity of 
the allowed range of the sterile neutrino DM mass, O(10 1 ) keV, which in our case is obtained 
from a very different approach. 

The Lagrangian of the right-handed neutrino sector, including gravity, reads (we use 
units h = c = 1 ): 

C = C-GR + £n r1 + Cy + £-1 (2.1) 

where 

CgR = 

Cn r i = V m N R i — -m N c riNri, 

Cy = 

Ci = -gvV,= -9vV„N R1 rN R1 , 

with R the Ricci scalar for the static spherically symmetric metric background 

5 V = diag(e y , -e\ -r 2 , -r 2 sin 2 6) , ( 2 . 6 ) 

where e u and e A depend only on the radial coordinate, r. The quantity m is the mass of the 
sterile neutrino, — !w(f[ 7 a 5 7 &] is the gravitational covariant derivative acting on a 

Majorana spinor, with the spin connection and [, ] the commutator. The right-handed 
sterile neutrinos Nr± satisfy the Majorana four-spinor condition, T c = T, together with 
'k = \k T C , where the conjugate spinor field T c = C'k and C is the unitary (C^ = C~ l ) 
charge conjugation operator, flipping the fermion chirality, i.e. (\k l) c = ('k c )/j is right-handed 
(R), whilst ('k L ) c = ('k c )j / is left-handed (L). The definition of chirality (handedness) is the 
standard one, 'S’l(R) = 7 ^lT 7 5 ) with the + (-) sign denoting Right-(Left)handed spinors, 

and 75 = Z 7 ° 7 1 7 2 7 3 , with 7 ^ the 4x4 Dirac matrices, satisfying 7 ^ 7 ^ + 7 ^ 7 ^ = 2 , where 
g^ is the (inverse) of the spherically symmetric space-time metric given above (( 2 . 6 )). 

The vector-meson mass is my, whose microscopic origin is not discussed here 4 , and 
Vfj, u = d^Vy — d u Vfj,, where the “Lorentz gauge condition” d ll V R = 0 has been applied for 
the vector-meson (VM) field V^. Notice that the massive-vector-mesons should not be 
viewed as gauge bosons if the fermions are Majorana. As is well known, the Lorentz gauge 

4 It may well come from an appropriate Higgs mechanism in the dark sector (with a Higgs field that is not 

necessarily the one of the SM sector). 


( 2 . 2 ) 

(2.3) 

(2.4) 

(2.5) 
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condition emerges then as a consequence of their equations of motion. Latin indexes denote 
flat tangent space indexes and are raised and lowered with the Minkowski r) a b metric. 

For simplicity we assume a minimal-coupling form of the vector field with the sterile 
neutrino current Jy in the interaction term Cj in the lagrangian density. This current 
is conserved if decays of sterile neutrinos are ignored. Such a coupling may also arise from 
linearisation of a Thirring-type four fermion vector current interaction Jy Jy^ by means of an 
auxiliary vector field A^ (which acquires dynamics upon implementing quantum corrections). 

In general one may add to (2.1) a Yukawa term, coupling the (three, in general) right- 
handed neutrinos to the active neutrino sector (see, e.g., the case of z^MSM [2-4]) 

£y u k = F a i£ a N RI (j) c + h.c. , 1 = 1,2,3 (2.7) 

where £ a are the lepton doublets of the SM, a = e, n, r, F a i are appropriate Yukawa couplings, 
and (j) c is the SM conjugate Higgs field, i.e. (f> c = zt 2 </>*, with T 2 the 2x2 Pauli matrix. 
Upon considering such a coupling, one obtains the stringent X-ray and BBN constraints of 
the mixing angle and mass of Nr ± depicted in figure 1, given that (2.7) implies decays of 
the heavy neutrinos Nj —>• vH, where H denotes the Higgs excitation field, defined via: 
4> = {(j>) + H. In such a case Jy is not conserved in time. However, in the context of ^rMSM, 
the lightest of the heavy neutrinos decay time is longer than the age of the universe, hence 
the latter can be considered as stable for all practical purposes, thuds playing the role of 
dark matter. 

For our purposes, as already mentioned, we concentrate here on this lightest neutrino 
and ignore such a mixing with the SM sector, setting F a \ = 0, in which case the lightest 
neutrino is absolutely stable. The important feature for us are the self-interactions of the 
right-handed neutrino, which will be used for ensuring phenomenologically correct values 
for the radius and mass of the galactic core. Since, as we shall see, the mass range we 
obtain is compatible with the one in figure 1, one may switch on the Yukawa term in a 
full phenomenological study, including the SM sector, and in particular neutrino oscillations 
and Early Universe physics (e.g. leptogenesis [2-4]), without affecting our conclusions. This 
stems from the very weak nature of the Yukawa couplings F a j as dictated by the seesaw 
mechanism which is assumed to be in operation here [6] so as to give a mass in the active 
neutrinos. For an order-of-magnitude estimate of such i^MSM (subleading) contributions to 
the effective four-fermion right-handed Majorana neutrino interaction strength we refer the 
reader to appendix C. 

A particularly interesting motivation to include coupling with the SM sector (active) 
neutrinos z/, is to be able to obtain a possible indirect detection method for the ‘inos’ through 
the decaying channel Nri —>• Z/+7, with a potential enhancement due to their self-interacting 
nature 5 . Particular attention should be paid to the recent observations by the Fermi satellite, 
providing evidence of a clear emission in the energy range 10-25 keV from the central region 
of the Galaxy [45]. The latter could find plausible explanation by means of a DM particle 
species with a mass of order 50 keV/c 2 , similar to the one obtained here. 

Notice that in eq. (2.4) we included a kinetic term for the VM-field. However, in the 
mean-field approximation we shall employ in this work, such kinetic terms are irrelevant, 
thus allowing contact four-fermion interactions among the right-handed neutrinos of Narnbu- 
Jona-Lasinio type to be studied in a similar way. In the latter case, the VM-field is auxiliary. 

5 In the context of the Yukawa term (2.7) such a decay pattern is obtained, e.g., from the decay of the 
Higgs to two photons. 
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From (2.1) one obtains the following equations of motion for the various fields: 

Gfj, v + = 0 , ( 2 - 8 ) 

+ m 2 v V v - g v J v y = 0 , (2.9) 

Nri + ^mN c R i = 0 , (2.10) 

where G^ v is the Einstein tensor and T^ v is the total energy-momentum tensor of the free- 
fields composed by two terms: Tj^ and Ty U , each of which satishes the perfect fluid pre¬ 
scription 

= (£ + V)u^u v - VgT , (2.11) 

with £ and V the energy-density and pressure which we define below. 

2.1 Relativistic mean-field approximation 

We now introduce the relativistic mean-held (RMF) approximation. In this approach, the 
system can be considered as corresponding to a static uniform matter distribution in its 
ground state 6 . Thus, the vector-meson field as well as the source currents are replaced by their 
mean values in this state, which, on account of space translational invariance, are independent 
of the spatial coordinates x; this and the requirement of rotational invariance imply that no 
spatial current exists, and only the temporal component of the current is non zero, i.e. 
Jy —>■ ( Jy ) = (Nri 7 °IV/ji) = (N R1 N R i). The last expression within brackets denotes the 
finite number density of right-handed neutrino matter times the temporal component of the 
pertinent (average) velocity. 

The RMF approximation allows one to solve the coupled system of differential equations 
(2.8-2.10) rather straightforwardly, to obtain directly the mean-field vector-meson as 

F 0 = -^r (2-12) 

m z y 

with the notation (Vo) = Vo and 

(Jv) = Jo =nu o> (2-13) 

where uq = e v G is the time-component of the (average) future-directed four velocity vector, 
and we have used the normalization condition u^u y = 1. 

The Major ana spinors in the RMF approximation can be simply expressed as the cor¬ 
responding momentum (Fourier) eigen-states with no x— dependent terms (see, e.g., [46]) 
\V(x) = 'F(fc) e ~ ik v xfl . Recalling that we are working here with a system comprising of a very 
large number N of fermions in thermodynamic equilibrium at finite temperature T, we can 
assume that the fermion number density is expressed in terms of the Fermi-Dirac distribution 
function f(k) 

n = e-G\N R1 (kh°N R1 (k)) = J d 3 kf(k). (2.14) 

where g is a spin-degeneracy factor for the Majorana spinors, the momentum integration is 
extended over all the momentum space, and f(k) = (exp[(e(fc) — fi)/(k R T)\ + l) -1 . Here 

6 As it is shown in section 3, this approximation is well justified when applied to all the fields (real and 
mediators) under the physical conditions of the quantum core, which is composed by a very large amount of 
fermions in a highly degenerate state, in some analogy with the physics of compact objects. 
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e(k ) = y/k 2 + m 2 — m is the particle kinetic energy, g is the chemical potential with the 
particle rest-energy subtracted off, T is the temperature of the heat bath, and kg is the 
Boltzmann constant. It is important to notice that we are working with the right-handed 
component of the full Majorana spinor d/, and so, although a full Majorana spinor (left plus 
right chiral states) is its own antiparticle implying a spin degeneracy g = 4, this is not the 
case for the singlet right-handed component Nri (viewed as a spin +1/2 fermion of one 
helicity state), for which g = 1. From now on we adopt this value for g. 


2.2 Thermodynamic equilibrium conditions and equations of motion 

We now introduce the thermodynamic equilibrium conditions. In the case of a self-gravitating 
system of semi-degenerate fermions at finite temperature in general relativity, in absence of 
any self-interactions (other than gravity) such conditions read [30]: e v ^ 2 T =constant and 
e u / 2 (g + m) =constant. The first equation corresponds to the Tolnran condition [47], and 
the second to the Klein condition [48]. In the presence of the vector-meson mediator in¬ 
teraction (2.5), it can be shown that only the Klein condition is modified; the generalized 
thermodynamic equilibrium conditions are (see, e.g., [49], for details) 

e v ! 2 T = const. , (2-15) 

e u / 2 (g + m) + gy Vq = e u ^ 2 (g + m + Cy n ) = const. (2-16) 


where the term gy Vo is interpreted as a potential energy associated to the new meson field 
V/. In deriving the middle equation of (2.16), we have used eqs. (2.12) and (2.13). 

We can then finally write the full system of Einstein equations (2.8) together with the 
thermodynamic equilibrium conditions (2.16) in the following dimensionless form' 


dM 

dr 

du 

dr 

dO 

dr 

f3 


47 Tr 2 £, 

^ M + Air'Pr 3 
f 2 ( 1 - 2M/f) ’ 


I i i Cym 2 - Lyra - q an 

1 du [ !+ 


Cymr n dn 


2/3 dr 
f3 0 e 2 


1 i Cvm 2 1 dh 

r -t- 47r 3 p d9 


(2.17) 

(2.18) 

(2.19) 

( 2 . 20 ) 


where the following dimensionless quantities were introduced: r = r/y, h = Gmy 2 , M = 
GM/x, £ = Gx 2 £, V = G\ 2 V, with m p = \J^/G the Planck mass, and we have introduced 
the dimensional factor x = 27r 3 ' /2 (l/m)(m p /m) with units of length, scaling as m~ 2 . We have 
also introduced the temperature and degeneracy parameters /3 = ksT/m , and 6 = ^/(fcgT), 
respectively; we have evaluated the constants of the equilibrium conditions of Tolnran and 
Klein at the center r = 0, which we indicate with a subscript ‘O’. We have also introduced 
the parameter 

Cy = gy/my, ( 2 - 21 ) 

which encodes information about the strength of the coupling of the effective interactions of 
the fermions (‘inos’) and the mass of the vector-meson mediator. The total energy-density 

7 For Cy = 0, the coupled system of differential equations (2.17-2.20) reduces to the standard form pre¬ 
sented in [30]. 
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and pressure £ and V contained in (2.11), can be split in two components, 


£ = £c + £v i V = Vc + Vy , 


( 2 . 22 ) 


with £c and Vc the contributions to the energy-density and pressure from fermions in the 
RMF approximation, calculated as {T 3 )n r1 = £c and (T 3 )n r1 = Vc respectively. They are 
fully determined by the distribution function f(k) (with particle helicity g = l) 8 


£c = m 


(2vr) 


V c = 

3 (2t r) 5 


m 

m 


i + 

i + 


e(fc) 

m 

e(fc) 

2m 


d 3 k , 
ed 3 k, 


(2.23) 

(2.24) 


while 

£v = Vy = -e _I/ myV r 0 2 = -Cyn 2 > (2.25) 

is the contribution from the VM-held. We shall next proceed to solve the system of equations 
(2.17-2.20), including a discussion on the boundary conditions appropriate for the description 
of the Milky Way, as a self-consistency check of the approach. 


3 Numerical solutions 


We now apply the theoretical formalism presented above to study the DM distribution on 
different astrophysical objects from spiral to large elliptical galaxies, for given boundary 
conditions in agreement with observations. At the end, we give the DM particle mass and 
total cross section constraints arising from the numerical analysis of the boundary-value 
problem. 


3.1 Spiral galaxies: the Milky Way 

The boundary conditions in this case are given by the request of the observational agreement 
of the inner quantum core and halo part with the following Milky Way properties: 1) the 
compactness of its ‘dark’ center (SgrA*), i.e. massive and compact enough to explain the 
dynamics of the S-cluster stars closest to the Milky Way’s galactic center, 2) the DM outer 
halo mass M^ and radius r^, and 3) the onset of flat galactic rotation curve with the specific 
value of the circular velocity Vh at r^. It is important to recall that we define the radius of the 
inner quantum core r c as the distance at which the rotation curve reaches its first maximum, 
and the outer halo radius 77 l at the onset of the flattening rotation curve, which occurs at the 
second maximum (see also figure 1 in Ref. [30]). Notice that the so called halo radius (and 
mass) represent the one-halo scale length (and mass) associated with the fermionic model 
here presented, and corresponding with the turn-over of the density profiles in total analogy 
as other halo-scale lengths used in the literature such as ro or r _2 as shown in figure 2. The 
rotation curve is given by the circular velocity 


v(r) 


GM{r) 


r - 2 GM{r)' 


(3.1) 


8 Alternatively, this contribution to the energy can be expressed as the expectation value of the energy 
(^yofco$), where E(k ) = fco are the energy eigenvalues of the corresponding Majorana Hamiltonian (see, e.g., 
[46]). 


- 10 - 













Following the above procedure, we shall constrain the physical conditions /3q and do, 
together with the physical parameters, such as the sterile neutrino mass m, as well as the 
coupling parameter Cy. We recall that the non-interacting case Cy = 0 of the model 
(2.1) has been recently solved in [29, 30], whose more general DM density profile shows the 
typical core-halo distribution composed of three different physical regimes as described in the 
introduction of the present article and demonstrated in figure 2. 



r (PC) 

Figure 2. Comparison of the Ruffini-Argiielles-Rueda (RAR) DM density profile [30] (non-interacting 
case Cy = 0) with the Navarro-Frenk-White (NFW) one [50], and with a cored Einasto profile [51, 52]. 
The RAR profile is here given for the parameters: (3q = 1.251 x 10~ 7 , 9q = 30 and m = 10.54 keV/c 2 . 
The NFW profile is PNFwif) = /9or 0 /[r(l + r/ro) 2 ] with p 0 = 5 x 10 _3 Mq pc -3 and r 0 = 25 kpc, and 
the Einasto profile is given by pe(t) = p - 2 exp [—2n(r/r_2) 1 ^ n — 1], with p_ 2 = 2.4 x 10~ 3 Mq pc -3 , 
r _2 = 16.8 kpc, and n = 3/2. The free parameters have been chosen to describe the typical properties 
of spiral galaxies [9, 10]. This picture was taken from the original version in [30]. 

Indeed, in figure 2 a solution with m ~ 10 keV/c 2 of the Cy = 0 non-interacting model of 
[30], is compared and contrasted with selected DM halo profiles used in the literature. In the 
sub-parsec core region and for core masses of r^j 10 6 M q typical of (Milky Way-like) galaxies, 
for an ino mass ~ 10 keV/c 2 , the thermal de-Broglie wavelength, A b = h/\f2TrmkBT, is 
larger than the inter-particle mean distance l of the inos, justifying the quantum-statistical 
nature of the core. A degenerate core with a very low temperature was found in [30] to be 
compatible with the outer halo properties such as the halo radius, mass, and rotation curves 
of order 10 2 km/s. In the Boltzmannian region, we have Xb/1 <C 1 and, as shown in [30], the 
specific value of the corresponding circular velocity in the flat region is intimately related to 
the temperature parameter, /3. 

We shall adopt here the ansatz that the self-interactions occur only in the quantum 
regime and thus within the core, where the thermal de-Broglie wavelength, 

Xb ~ y/2irmk B T' ^ 

is larger than the inter-particle mean distance l, i.e. Xb/1 > 1- The reader should recall that 
at this quantum regime, two-particle interactions were predicted to be non-negligible [23]. 
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To this end, we set: 


C v (r) 



at r < r m when Xb/1 > 1, 
at r ^ r m when Xb/1 < 1, 


(3.3) 


where Co is a positive constant and r m = r c + Sr is the core-halo matching point, with r c 
the core radius and Sr the thickness of the core-halo intermediate layer. As we shall show, 
Sr -C r c , and thus the core-halo matching satisfies r m ~ r c . In the regime r ^ r m , where the 
DM distribution is in a much more dilute state (i.e. Xb/1 1), there is the transition from 
the quantum degenerate state to the Boltzmannian one. 

As we show below, the density profile obtained for the interacting case has a similar 
behavior, with the aforementioned three different regions, as the non-interactive case Cy = 0 
[30]. We normalize hereafter the coupling constant Co, for the sake of reference, to the 
Fermi constant of the SM weak-interaction, i.e. we introduce the dimensionless constant 
Co = Cq/Gf■ We define the SM Fermi constant only for normalization purposes, thus Co 
must not be thought as a fundamental interaction strength (i.e weak) of the SM. Indeed, 
the fact that the effective interactions considered here are mediated by a chargeless VM 
field playing the role of neutral-current interactions through the scattering channel, implies 
that the inos remain unaffected except for momentum transfer. Therefore, we here adopt 
a complete phenomenological analysis by studying the maximum possible range of effective 
interactions strengths which are in agreement with the Milky Way observables. 

The mass M c of the degenerate quantum core must agree with the mass enclosed within 
the region bounded by the pericenter of the S2 star. At the same time, we use the pericenter 
of S2 as an upper limit to the core radius r c ( S2 ), i.e. [36, 37] 

M c = 4.4 x 10 6 M 0 , r c(S 2 ) = 6 x 10~ 4 pc. (3.4) 

There is an error of 8% in the above value of M c due to the uncertainties in the measurement 
of the distance to the galactic center Rq = 8.33 ± 0.35 kpc, while the error in the pericenter 
of the S2 star is of about 4% [37]. The above parameters imply a central density of order 
~ 10 le) Mo/pc 3 , which is almost five orders of magnitude larger than the one obtained for the 
model without self-interactions [30] with the same core mass. It is important to make clear 
that any core radius rs c h < r c < r c ( s2 ) is accepted within our phenomenological treatment, 
implying central densities in the range 10 16 M Q /pc 3 < po 10 23 Mq/pc 3 , with rs c h the 
Schwarschild radius of a black hole of 4.4 x 10 6 Mq. Indeed, as we show below, already for an 
ino mass m ~ 350 keV/c 2 it is possible to obtain a critical core DM core of fully degenerate 
inos of mass M c = 4.4 x 10 6 Mq with a radius r c ~ 2.5rg c h. The critical objects are the last 
equilibrium configurations, just before undergoing gravitational collapse (see also Ref. [53]). 

For the observables in the halo region we adopt the fitting procedure outlined in 
ref. [54]. According to that work, the DM best-fit distribution for the Milky Way is pro¬ 
vided by the two-parameter Burkert profile with a specific central density parameter p° B = 
2 x 10 _2 Mq/pc 3 , and a dark halo length scale parameter h = 10 kpc. For our fermionic 
model, this corresponds to a halo radius r^, defined at the maximum of the corresponding 
rotation curve at the onset of the flat behaviour, and leads to an associated halo velocity 
Vh and mass Mf l given by (the reader is invited to observe the excellent matching between 
Burkert and fermionic profiles around in figure 3): 

rh = 32.4 kpc , Vh = 155 km/s , Mh = 1.75 x 10 11 M Q , (3-5) 
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where the subscript h indicates quantities at the halo radius. All the halo parameters are 
subject to an error of ~ 10% [54], The above value of the circular velocity determinates 
the value of the temperature parameter at the halo, /3 q. For these parameters, we obtain 
$ = 1.065 x 10" 7 . 

We discuss now the core-halo transition. There, the generalized Tolman and Klein equi¬ 
librium conditions have to be fulfilled. The Tolman’s condition together with the condition 
imposed by the continuity of the spacetime metric, lead to the continuity of the temperature 
parameter /3(r m ) = /3 q. Now, from the Klein’s condition we can obtain the jump in the 
degeneracy parameter at the matching point r m , where the (diluted) halo region begins: 



n(r m ) = ^ m3 ff )3/2 (T 1/2 + $F 3/2 ) , (3.7) 

where 6 g is the value of 9 from the halo side, and the generalized Fermi-Dirac integrals are 
evaluated at r m : Fj = dxx J (1 + Pqx/2) 1 ! 2 /[I + e x ~ 6l ( rm )]. 

3.2 Other galactic structures: large elliptical galaxies 

In analogy with the Milky Way case, we now apply our SIDM model to other galactic 
structures, such as large spiral and elliptical galaxies, where clear evidence exists for massive 
BH-like structures at their centre, together with a DM halo counterpart. 

In order to give a universal explanation for the galactic DM in terms of an ino mass 
and an interaction constant (or cross-section), we next apply our theory to larger galaxies for 
m = 47 keV, and give the possible values of Cy (2.21) in agreement with the different galaxy 
observables. We proceed with two different kinds of typical elliptical galaxies each harboring 
a different characteristic dark massive object at the center: 

(i) M c = 2.3 x 10 8 Mq Elliptical, 

(ii) M c = 1.8 x 10 9 Mq Large Elliptical, (3-8) 

both contained within sub-pc scales [55]. Notice that the above cases are representative 
examples and other intermediate cases between normal spiral galaxies and the large elliptical 
ones are also contained among the possible solutions of our model (see section 3.3), but are 
not given explicitly here for the sake of brevity. 9 * 

The halo parameters are chosen from the observationally inferred correlation between 
central mass concentrations and dark halo masses ( M c -Mh ) as obtained in [56] (see figure 
5 there). This is in analogy with the study of [30], but in our SIDM case we do reach the 
upper end of the correlation. Thus, typically, we should have 

(i) Mh = 4.1 x 10 12 Mq , rh = 60kpc, Elliptical, (3-9) 

( ii ) Mh = 1.1 x 10 13 Mq , r/ l = 85kpc Large Elliptical. (3.10) 

The values given above were considered at the one-halo scale-length of our model (located 
at the maximum of the halo rotation curve), which is similar to the NFW halo scale-length r s 

9 Regarding the applicability of our approach to dwarf galaxies, we remark that the rather low central 

degeneracy values in the phase-space distribution of these systems [30] lead to not very massive cores, therefore 

not exhibiting massive BH-like features (i.e. with masses M c < 10 4 Mq for m ~ 10 1 keV). Hence, they do not 

seem to constitute interesting cases for a detailed study of the effects of self-interacting DM. 
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Milky Way (M c = 4.4 x lO a M 0 ) 


m (keV) 

Co 

o 0 

P 0 

r c (pc) 

dr (pc) 

0{r m ) 

47 

2 

3.70 x 10 3 

1.065 x 10" 7 

6.2 x 10” 4 

2.1 x 10" 4 

-29.3 


10 14 

3.63 x 10 3 

1.065 x 10" 7 

6.2 x 10" 4 

2.2 x 10" 4 

-29.3 


10 16 

2.8 x 10 3 

1.065 x 10" 7 

6.3 x 10" 4 

2.4 x 10" 4 

-29.3 

350 

1 

2.40 x 10 6 

1.431 x 10" 7 

1.3 x 10” 6 

6.7 x 10" 7 

-37.3 


10 14 

1.27 x 10 5 

1.104 x 10" 7 

5.9 x 10" 6 

9.4 x 10" 7 

-37.3 


4.5 x 10 18 

1.7 x 10 1 

1.065 x 10" 7 

5.9 x 10" 4 

2.0 x 10" 4 

-37.3 



Elliptical (M^ r = 2.3 x It 




47 

2 

1.76 x 10 5 d) 

1.7 x 10" 6 

7.9 x 10" 5 

3.9 x 10" 5 

-31.8 


10 14 

5.8 x 10 4 

1.4 x 10" 6 

1.4 x 10" 4 

4.8 x 10" 5 

-31.8 


10 16 

1.5 x 10 4 

1.3 x 10" 6 

3.0 x 10" 4 

7.0 x 10" 5 

-31.8 

Large Elliptical (M c = 1.8 x 10 9 Mq) 

47 

10 16 

1.02 x 10 4 

3.0 x 10" e 

3.8 x 10 -4 

1.8 x 10" b 

-32.8 


Table 1. Set of model parameters for three different galaxy types analyzed that solve the correspond¬ 
ing boundary-value problem imposed by the given galaxy observables, from spiral (Milky Way) to 
large elliptical galaxies, d) Critical central degeneracy parameters (0g r ), associated with the turning 
point or last stable equilibrium solution [53]. 


(r h ss r s /0.6), where the halo mass was originally obtained in [56]. These halo magnitudes, 
combined with the DM halo morphology of our model, imply typical halo velocities of Vh = 
540 km/s in case (i), and Vh = 730 krn/s in case (ii). With these values, we finally obtain 
(analogously as done for the Milky Way in appendix A) the following temperature parameters 
at the halo: /3ft = 1.3 x 10 -6 , /3ft = 2.4 x 10~ 6 for (i) and (ii) respectively. The parameters 
of the model to be used for solving the above boundary conditions, are as follows: we keep 
the ansatz (3.3) for the interaction constant, while we set m = 47 keV. 

3.3 Novel DM mass constraints 

Following the above procedure, we summarize in table 1 the solution of the boundary-value 
problem which fulfills the core and halo observables from Milky Way (3.4), (3. 5), and Elliptical 
and large elliptical galaxies (3.8)-(3.10) respectively. The calculations were done for the 
maximum allowed possible range of the interaction constant Co, central degeneracy 0q and 
ino mass m. Even if the upper limit in the sterile neutrino mass (m < 50 keV/c 2 ) is imposed 
by cosmological and astrophysical constraints under the assumption of mixing with the SM 
sector (c/. figure 1), we also explore larger (phenomenologically) values of the ino mass, which 
is possible for sterile neutrinos that do not interact or have negligible interactions through a 
Higgs portal with the active sector. 

Two important conclusions can be drawn from the numerical analysis presented in 
Table 1: 

I) For m < 47 keV/c 2 and m > 350 keV/c 2 there is no pair of parameters (Co,0o) able 
to be in agreement with the Milky Way observables. While m = 47 keV/c 2 is the lowest 
admissible particle mass up to which the core observational constraints are fulfilled (within 
observational errors), m = 350 keV/c 2 is the uppermost bound set by the reaching of the 
critical core mass for gravitational collapse [53], M£ r oc /to 2 ~ 4.4 x 10 6 Mq, where M p i is 
the Planck mass. For m = 47 keV, one reaches the critical mass M£ r ~ 10 8 M@ corresponding 
to a massive BH alternative for elliptical galaxies (when Vy < Vc)\ but it is also possible to 
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reach core mass values as large as M c 10 9 Mq (when V\/ r^i Vc), and applicable to massive 
BH alternatives in large elliptical galaxies. In the latter case the quantum core can reach 
radii as small as r c ~ 2rg c h (see last line in table 1), developing interior sound-wave speeds 
as large as 5% of the speed of light. 

It is worth to notice that for all the galaxy types analyzed, there is a common window 
of interaction coupling constant parameters (Co), the value Co = 10 16 constitutes a very 
interesting case, allowing for a successful universal application of the model all the way from 
spiral to large elliptical galaxies. 10 

II) As the value of the coupling constant Co increases from unity, the contribution to 
the total energy and pressure from the meson-vector field (~ Con 2 ) becomes more and more 
relevant. For instance, as can be seen in table 1 in the Milky Way case, for Co ~ 10 14 
and for m = 47 keV, a slightly lower value for the central degeneracy is needed to have the 
same core mass as compared with the Co ~ 1 regime. In other words, if the same central 
degeneracy as in the former Co ~ 1 case is used, an increase of ~ few % in the core mass M c 
would appear. For this lower ino mass bound, the self-interactions cannot exceed Co ~ 10 16 , 
because otherwise the now lower central degeneracy needed to compensate for the core mass, 
would be too low to fulfill with the upper core radius constraint r c ^ 2 )- More evident is 
the case when the ino mass reaches m = 350 keV/c 2 , where the highest interaction regime 
Co rsj 10 18 fulfilling the core radius and mass, is reached at a central degeneracy about two 
orders of magnitude lower with respect to the Co = 1 case. 

3.4 Cross-section constraints 

It is possible to establish, within our theoretical approach, a direct link between the total 
cross-section a and the interaction strength Co = ( gy/my ) 2 - This will allow us to compare 
our results with the ones given in the literature, regarding the total cross-section per DM 
mass, a/m. To this end, we shall consider a four-fermion (elastic) scattering for the inos, 
with a massive-vector boson V^ as mediator field. 

The total cross-section in the center-of-mass (CoM) system for two incident right-handed 
Majorana neutrinos Nr \ with four-momentum (pi = (E,p/)), i = 1,2, (with E 2 = p 2 + m 2 ) 
which collide through an elastic scattering picture, and produce two final identical particles 
with momentum (p' t = (E 1 , p( = Pi)), i = 1,2 (E’ = E) is given by (see appendix B for 
details) 

*CoM = (—) 73^[29 m 2 + 89 p 2 + 89/3^^] . (3.11) 

\my / 4^7r z p z + m z 

where we have used 9' w = 0 since we only have one massive-vector mediator. 

We now calculate the total Nri-Nri scattering cross-section in the quantum core of 
the Galaxy. For this we use the following approximations leading to a simplified version of 
eqn. (3.11). In a typical quantum-core (non-relativistic) one-particle momentum p is given 
the Fermi momentum (at the core of the configuration) p ~ pp = (3ir 2 h 3 pc/m) 1 / 3 . For typical 
core densities used in this work, p c ~ 10 l ®~ 23 Mq /pc 3 , and m = 47-350 keV, which leads 
directly to the following ‘low-energy limit’ p 2 <C m 2 for our particles. This is the opposite 
with what one generally finds in laboratory collision-particle experiments (i.e a ‘high-energy 
limit’), which is easily understood because the sterile neutrinos are in a very low temperature, 

10 If, instead, the interaction constant is forced to agree with the N-body simulation results for the total 
DM cross-section (i.e. Co = 7 x 10 8 for m = 47 keV, as detailed in section 3.4), then the applicability of our 
model is reduced up to elliptical galaxies with dark compact cores of ~ 2 x 10 8 Mq. 
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and a high degenerate regime in the core. With all this, equation (B.23) reads 

core ~ { 9 V {Z Vr ^ 2 ( P 2 /m 2 « 1) ■ (3.12) 

Equation (3.12) links the (dimensionless) interaction constant of our inos (expressed relative 
to the (weak interactions) Fermi constant Gf = 1.166 x 10~ 5 GeV~ 2 ), (cf. (2.21)), 

C v =(—)g- f \ (3.13) 

\m v J 

with the total cross-section and the particle mass. Thus, if we constrain the total cross- 
section to the N-body simulation value a tot /m = 0.1 cm 2 /g [43], our coupling constant Cy 
would be constrained to the value 

Cy € (2.6 x 10 8 ,7 x 10 s ), (3.14) 

for ino masses in the range m € (47,350) keV. It worths noticing that for Cy 10 s Gf, the 
mass of the massive-vector meson would be constrained to values my < 3 x 10 4 keV, in order 
to satisfy gy < 1 as requested by the self-consistency of the perturbation scheme we have 
applied to compute the cross-section. 

We can further try to get an absolute lower bound for Cy. Interestingly, it can be 
obtained by answering the question as to which physical conditions need to be fulfilled by 
the keV particles in the galaxy for a self-interacting DM regime to appear. A conservative 
answer one might give is that a should be large enough so that a scattering probability among 
the inos should occur at least once during the age of the galaxy (t age ), that is, the product 
of the scattering-rate per particle (T) times t age be larger than unity: T t age > 1. For this 
we first calculate T, which is linked to the cross-section a via [57] 

T = cr\v re i\n, (3.15) 

where v re i is the relative velocity of the interacting inos and n the particle number density. 
The above formula can be written as an order of magnitude expression as follows: T ~ 
a/m 2 p po, with p a typical momentum of the inos in the quantum core (i.e. p = pp) and 
Pq the central density (valid in the low energy regime approximation). Then, by assuming 
typically t age ~ 10 16 s (i.e. redshift 2 ~ 10 at galaxy formation epochs), one obtains, for 
m = 47 keV and po 10 16 M Q /pc 3 (as for the Milky Way case, see table 1): 

a/m > l(T 18 cm 2 /g, (3.16) 

directly implying from our cross-section formula within the low energy approximation (3.12), 
that Cy > 2 Gf, that is the interaction strength can never be smaller than the weak interac¬ 
tion Fermi coupling. 

4 Discussion 

It is interesting to notice that the degenerate ~ 10 1 keV fermion core can reach core radii 
small enough to be suitable for the SgrA* observational constraints, as well as to reach 
BH-like compactness as in the case of larger elliptical galaxies. This alternative approach 
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Figure 3. Left: mass density profiles for to = 47 keV/c 2 in the interaction regime Co = 10 16 
where core and halo Milky Way observational constraints (3.4-3.5) are fulfilled, compared with the 
non-interacting case (Co = 0) for the same ino mass in disagreement with the core observables. 
We also show for comparison the two parametric Burkert profile pb/[{ 1 + f/h)( 1 + (r/h) 2 )] with 
Pb = 2 x 10 _2 Mg/pc 3 and h =10 kpc, which is the best DM halo fit of the Milky Way according to 
[54]. Right: degeneracy parameter profile in the interaction regime Co = 10 16 for the same ino mass 
as in the Left panel. 
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Figure 4. Left: mass density profiles for m = 47 keV/c 2 in the same interaction regime as in 
figure 3, where the large elliptical core and halo observational constraints (3.8-3.10) are fulfilled. 
We also show for comparison the two-parameter Logaritmic profile pg[l + l/3(r/r(i) 2 ]/[l + (r/r'f) 2 } 2 
with P q = 2.0 x 10 ~ 2 M Q /pc 3 and = 35 kpc (r^ ft! 2.4r£), usually used in typical large elliptical 
galaxies [58]. Right: degeneracy parameter profile in the interaction regime Co = 10 16 for the same 
ino mass as in the left panel. Notice that the small bump in 6(r) around r = r c is absence in the 
non-interacting case, and thus it is originated from the presence of the meson-field in the system of 
equations (2.17-2.20) 


acquires special interest for ongoing and future observational campaigns (e.g. the BlackHole- 
Cam project 11 ), which would allow to verify the general relativistic effects expected in the 
surroundings of the central compact source in SgrA*; leading to a deeper scrutiny for the 
not-yet confirmed black hole hypothesis. 

In figures 3 and 4 we show, for comparison, the overall density distribution p(r) cor- 

11 http: //horizon-magazine, eu / space 
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responding to the Milky Way and a large elliptical galaxy respectively, for a common self¬ 
interaction constant. For the Milky Way we show as well the p(r) profile in the non-interacting 
case for the same ino mass m = 47 keV/c 2 . This comparison shows that, while in the non¬ 
interacting case (Co = 0) the core observables (3.4) are not fulfilled, the presence of self¬ 
interactions allows to have higher degenerate cores satisfying both the core and halo Milky 
Way observables (3.4-3.5). It is important to notice that the density profile in the observa- 
tionally well constrained halo region of figure 3, coincides with the one obtained in ref. [30] 
in the absence of self-interactions for the same ino mass, and are in good agreement with the 
Burkert profile, best DM halo fit for the Milky Way, as shown in [54]. 

At this juncture we should point out that in our analysis above we did not discuss explic¬ 
itly the role of baryonic matter, which of course is mainly concentrated through the central 
bulge and disk regions of galaxies. Its inclusion does not change the important conclusions 
of our work that the introduction of WDM fermion self interactions affects the core/halo 
structure and in particular induces higher central degeneracies and higher compactness of 
the inner quantum core. The key result presented here as well as in [30] is that the DM 
contribution is predominant in the inner core (below sub-pc scales), and in the halo region 
at the onset of the flat part of the given rotation curve; while in between baryonic matter 
prevails. Indeed, we can see from figure 3 (left panel) that for the Milky Way, our model cor¬ 
rectly predicts both the value and flattening of the circular velocity at distances r > 10 kpc. 
A more complete quantitative analysis, including baryonic matter, is left for a future work. 

We would also like to make one last, but not least, observation regarding the range of 
the self-interacting ino masses, m > 47 keV/c 2 . If we identify the inos with the (lightest) 
right-handed neutrino of the z^MSM model [2-4], then the latter must have a very weak 
mixing angle with the SM lepton sector, and its mass must be less than 50 keV/c 2 , otherwise 
the model would not be consistent with the current phenomenology, as can be seen from 
figure 1. The above considerations, then, leave a very narrow range of the self-interacting 
‘ino’ mass 47 < m < 50 keV/c 2 , for the right-handed neutrino to play both a role as a WDM 
candidate and a provider of a core-halo galactic structure in accordance to observation. Such 
constraints are of course alleviated if any mixing of the ino with the standard model sector is 
avoided, as done in the current article. Nevertheless, we find quite intriguing the fact that, 
starting from two entirely different approaches, one from particle physics, and the other from 
pure galactic astrophysics, one finds a consistent regime of ‘ino’ masses within the WDM 
range. We believe that this is not a coincidence, and the aim of the current paper was to 
alert readers from these different communities to this important fact. 

Finally, before closing, it is appropriate, as announced in the introduction of the article, 
to place our fermionic keV DM approach in context with the current state of affairs of 
cosmological DM and structure formation and some of the important open issues still faced 
by the ACDM cosmology [59], such as: (i) the core-cusp problem [41], (ii) the “lost satellite” 
problem [42] and (iii) the, so-called, “too big to fail” problem [60, 61], the latter being a 
discrepancy between the most massive subhaloes arising within CDM and the dynamics of 
the brightest dSph galaxies of the Milky Way. All these problems have their root in the 
fact that cold DM particles have too short free streaming length during the epochs of galaxy 
formation, and therefore they form too clumped and too many structures than those observed. 

Our model provides a natural solution for (i) because the density profiles based on 
fermionic phase-space distributions develop always an extended plateau on halo scales (start¬ 
ing after the quantum core), in a way that resemble Burkert or cored Einasto profiles [30] (see 
figures 2 and 3). Regarding the issues (ii) and (iii), it is important to bear in mind that our 
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model does not directly deal with structure formation mechanisms, nor has employed (as yet) 
numerical simulations at such scales, and hence we are not in a position to make any concrete 
statement on these two issues. Nevertheless, for particle masses in the few 10 1-2 keV range 
as obtained here, it has been extensively shown by now that the behaviour in the power 
spectrum (up to ~ Mpc scales) is practically indistinguishable from that of standard CDM 
particles [4], thus maintaining the expected results from large-scale structure observations. 
There are issues, such as re-ionization, that we have still not examined, and thus at this stage 
we cannot make any concrete statements as to how much our model provides a substitute for 
CDM, although several of its features, as we have explained above, are indistinguishable from 
it. In fact, it may well be that there exist more than one DM species in the universe, and in 
this respect, our self interacting warm (right-handed neutr)“inos” play an important role in 
the galactic-core structure, which was analyzed above, but a complete explanation/resolution 
of the large-scale structure problem in the cosmos may require synergies among different DM 
species, including CDM and WDM. More work is needed to arrive at firm conclusions in 
these matters. 

In this respect we mention for completeness that several proposals have been made re¬ 
cently towards a resolution of these issues within standard ACDM N-body simulations, includ¬ 
ing the self-interacting DM approach of [43], as well as baryonic feedback processes [62, 63]. 
Moreover, within the realm of (Newtonian) N-body simulations, AWDM cosmologies based 
on particles of few keV have been developed to tackle the aforementioned ACDM discrepan¬ 
cies (i)-(iii) [64, 65]. Particles with such a small mass can suppress structure formation on 
small scales, due to their larger free streaming length caused by appreciable thermal veloci¬ 
ties. Nevertheless, the fact that recent observations suggest that the number of Milky Way 
satellites is an order of magnitude greater than that predicted by WDM numerical simula¬ 
tions, casts doubts for the few keV WDM scenario [66]. Moreover, such light (1 — 3 keV) 
particles are also in strong tension with actual lower keV bounds set by current Ly-a forest 
constraints [6, 67]. 

All these issues present serious challenges for N-body simulation-based cosmologies (i.e 
AWDM in this case), associated with the ‘too warm’ nature of the particles involved. In¬ 
terestingly, the fact that the particle mass in our model is ‘colder’ by a few keV compared 
to those WDM models, implies that our model does not suffer from such standard WDM 
problems. This, together with the fact that the model tackled successfully the important 
core-cusp discrepancy, as mentioned above, and that it avoids several of the undesired halo 
features characterizing the ACDM paradigm, offers significant support to the idea that our 
self-interacting model of right-handed neutrinos provides a physically important DM species, 
which may co-exist harmonically with other DM structures in the universe. In this respect, we 
believe that complementary searches for such keV right-handed neutrinos, either in neutrino 
oscillation experiments or elsewhere, are important. 
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A Central temperature parameters 

A sufficiently precise determination of the central temperature parameter /3 q in the low 
relativistic regime of the model, when applied to normal galaxies, can be understood through 
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the following two concepts (the value of the speed of light c is here given in krn/s): 

1) Boltzmann regime at r ~ r m : Since at r > r m the degeneracy parameter fulfills 
9(r) -C — 1, the Fermi-Dirac statistics necessarily approaches the pure Boltzmann regime. 
The Boltzmann distribution function is characterized by the familiar one-dimensional velocity 
dispersion, a, which is independent of the radius 

a 2 = kBT/m. (A.l) 


2) Classical isothermal-sphere condition. A classical self-gravitating system of Boltzmann- 
like particles in hydrostatic equilibrium is described by the isothermal-sphere model. The 
relation between the circular velocity v c (r) and a for an isothermal-sphere model is v 2 (r) = 
—a 2 (d\np(r)/d\nr), where p(r ) is the mass density (see, e.g., [7]). Different cored solutions 
to p(r) depend only on the constant initial conditions p$ and a, implying a universal behav¬ 
ior (scaling) of the density profile. Thus, the logarithmic derivative evaluated at the halo 
radius rh (defined at the maximum of the velocity curve, i.e. the onset of the flat part) 
is (din p(r) / d\nr)\ Th = —2.51. This implies v\ = 2.51c 2 , and, hence, using eq. (A.l), one 
obtains 


08 



(A.2) 


which for Vh = 155 km/s gives 


^ = /3(r m ) = 1.065 x 1(T 7 . 


(A.3) 


This is the value we use in our phenomenological analysis of section 3. 

Finally, notice that, in order to obtain the central temperature parameters j3o appearing 
in table 1 , we use the Tolman condition for isothermality e u / 2 T = constant. The latter, 
together with the definition of /3 = ksT / (me 2 ), implies the relation (2.20), which expresses 
the temperature parameter at any given radius in terms of the central temperature parameter 
/3o- For example, in case the massive quantum core of SgrA* has a small compactness (i.e. 
r c ~6x 10~ 4 pc), as dictated by GM c /r c ~ 10~ 4 , one necessarily has the following condition 

UQ-u(rm) 

for the metric factor between core and halo e 2 ~ 1. Consequently, the following 

relation can be established f5 q = /3o at three-digit precision. Instead, in cases with higher 
core compactness (i.e. when r c ~ 10~ 6 pc in the Milky Way case, or when r c ~ 10 -5 pc in 
the elliptical galaxy case) as shown in table 1, slightly higher values of (3 0 are obtained by 
the use of the Tolman condition. 


B Total cross-section a tot within a four-Majorana-fermion (elastic scatter¬ 
ing) interaction 

A scattering process consisting in two incident particles with four-momentum (pi = (Ei, Pi)), 
i = 1,2, which collide and produce two final particles with momentum (p( = (E(, pQ), 
i = 1,2 (all with definite polarization/spin states) has a total cross-section in the center-of- 
mass (CoM) system given by (see eqn. (8.19) in [68]) 

««-rw 2 ’-5/"- (Bi> 

The integration in the last member is made over the complete solid angle, while the 
integrals involving the (0,<j>) scattering angles were performed over the forward hemisphere 
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Figure 5. Feynman diagrams for v s - v s scattering, with Vq the massive-vector boson mediator 


(0 < 0 < tt/2) corresponding to physically distinguishable events, and proper for a process 
with two identical particles in the final state, as it is our case of interest. With ( da/dQ)c 0 M 
the differential cross-section in the Center of Mass system given by the general formula [68] . 
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where the sub-index 1 runs over all external leptons in the scattering process, and M is the 
Feynman amplitude of the process. 

A lepton-lepton (y s - v s in this case) scattering process has two possible Feynman 
diagrams (differing by the interchange of the final lepton states 1' -f -1 2') given in figure 
1, with two different associated Feynman amplitudes M a and Mb, built up below by the 
use of the Feynman rules associated to each diagram (a) and (b) respectively (see e.g. [68], 
appendix B). It is important to realize that being the sterile neutrinos neutral (no electric 
charge), the only possibility for the boson mediator in the interaction (at lowest order in 
the interaction term of the Lagrangian), is to be a neutral vector boson Vo. This reasoning 
is implemented in complete analogy to the scattering of leptons within electroweak theory, 
where the Zq boson can be viewed as the analogue of our dark-sector Vo boson. 

The Feynman amplitudes M a and Mb read, in accordance with the corresponding 
Feynman rules, as follows 


2 

M a = ~ 16cos^) [^i(Pi)7 a (l - ^)u ri {jPi)]iD Fap (k = p 2 -p / 2 ,m)[u S2 (p / 2 )7 a (l - 7sK 2 (P2)] 
Mb = 16co S ^) [^2(P2)'V a ( 1 “ l^) U rAVl)\iD Fa p{k = P 2 -p'i,m)[u si (p'i)7 a (l - 75K 2 (P2)] 

Where gy is the Yukawa coupling constant of our (only) interaction term in the la¬ 
grangian (£/), with 6' w the ‘pseudo-weak’ mixing angle 12 , u r ( p) (ii s (p’)) the four-Majorana 
spinor corresponding to a particle of momentum p, positive energy E, and spin indexes 

12 This can be formally implemented within our (minimally extended) beyond SM sterile neutrino, through 
the linear combinations (dark sector fields) of B^x) (e.g. [68] Eqtn. (12.45)), in total analogy with the 
electroweak mixing angle 9w- 
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r = 1,2 (s = 1,2) of the neutrinos in the initial (and final) states respectively (being u the 
adjoint). And Dp a ^{k,m) the massive-vector boson propagator which reads 


D F a p(k, m) 


-g a p + kgkp/m\ 

k 2 — rriy + ie 


(B.3) 


with g a p the metric tensor with signature mostly minus and locally Minkowskian (equiv¬ 
alence principle), my the vector boson mass and e a very small parameter (to avoid pole 
divergences). As it will be clear later, the typical sterile neutrino (one particle) energies 
we are dealing in our galactic environment, implies we are in a ‘low energy limit’, and we 
can safely assume k 2 <C my (for ‘normally’ heavy bosons Vo with masses above keV, as we 
assume here). With this inequality we have for (B.3) 


D F (k,m) ^ {k 2 /m 2 v 0), (B.4) 

my 

which leads to the following simplifications in the Feynman amplitudes A4 a and Xl b 

M a = (^) lecos^ tMPib^ 1 - 75)«r 1 (pi)]['w S2 (p / 2 )7a(l - 75K 2 (P2)] 

Mb = i 16cos^(( 4) [^(Pab 0 ^ 1 - 75K 1 (pi)]fosi(p / i)7a(l - 75K 2 (P2)] 

The Feynman amplitude factor in (B.2) can be written now in terms of the above 
amplitudes \M \ 2 = \M a + M b \ 2 = \M a \ 2 + |-M fe | 2 + M a M* b + (M a M b )*. With * the 
complex conjugate. 

In following, we assume that the ‘polarization’ state of our (initial and final) neutrinos 
is not known, and therefore we shall calculate the (standard) ‘unpolarized’ cross-section, for 
which we must average \JV[\ 2 over all initial and final lepton spins. Indeed, notice that our 
Majorana singlets Np\ are of right-handed chirality, but nevertheless they may have both 
positive and negative helicities. This averaging procedure consists in average over initial 
spin states (i.e. 1/4^) rl ^/] r2 ), an< ^ summing over final spin states (]T) sl ( see e -S- [68], 

section 8.2). Therefore, each Feynman amplitude term (\M . a \ 2 ,\M . b \ 2 ,etc.) is now replaced 
by 

X = \ E E E E I M(rl,r2, si, s2)| 2 , (B.5) 

rl r2 si s2 

leading to a differential cross-section (B.2) of the form 

(^) <x (X aa + X bb + X ab + X* ab ), (B.6) 

\auj CoM 

where 


X aa 

— I v \M I 2 

4 spins 1 1 &I 

(B.7) 

x bb 

= 3 Espins \Mb \ 2 

(B.8) 

x ab 

= 3 EspinsMaMl 

(B.9) 
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By inserting equations (B.5) in the above Feynman amplitude factors, we have corre¬ 
spondingly 


Ki(Pib a (l -75)«n(Pi)] 2 EE Ka(P2)7a(l - 75K 2 (P2)] 2 

rl si r2 s2 

x bb = c EE K 2 (P2)7 a (l -75Ki(Pl)] 2 EE K 1 (Pl)7a( 1 - 75)«r 2 (P2)] 2 

rl s2 r2 si 

x ab = -c EEEEi* 51 (p'l)7“(l - 75 Ki (pi)] [u S2 (p' 2 )7a(l - 75)«ra (P 2 )]® 

rl r2 si s2 

[«5 2 (P2)7“(1 - 76)“r 1 (Pl)][“* 1 (p , l)7«(l - 75)«ra(P2)] > 

with C = (gy/my ) 4 1 /[4 5 cos 4 (0' w )\. Then, by using the Hermiticity condition of the 
gamma matrices (( 7 ^)^ = 7 0 7 ^ 7 °) and the regrouping spinor-component technics (e.g. [ 68 ], 
section 8.2), it can be shown that the Feynman amplitude factors written above can be easily 
expressed in terms of traces of products of gamma matrices as follows: 


X aa = CTr[A+ (pib“(l - 7 5 )A+(pi) 7 /3 (l - 7 5 )] Tr[A+(p' 2 ) 7a (l - 7 5 )A+(p 2 ) 7/J (l - 7 5 )] 
X bb = <7Ti-[A+(p' 2 ) 7 “(l - 7 5 )A+(pi) 7 /3 (l - 7 5 )] Tr[A+ (pi)7 a (l - y 5 )A+(p 2 ) 7 /3 (l - 7 5 )] 
X ab = -C"Tr[A+ (pi)7 a (l - 7 5 )Am(Pi)7 /3 (1 ~ 7 5 )Am(P 2 ) 7«(1 ~ 7 5 )Am(P2)7 /?( 1 ~ 7 5 )] 

where A^(p) = A+j(p) = Yll=i u ra(p)u r p(p) = is the Majorana (positive) 

energy projection operator (see next section for details and difference w.r.t the Dirac pro¬ 
jector), with p = 7 A pa -In next we rewrite the above traces as X aa = (C/m 4 )A a/3 B a p and 
X bb = (C/m 4 )E a P F a p and proceed to calculate each factor-element A, B, E, F. 

A *P = TrK-p^ + m/ 2) 7 q (1 - 7 5 )(-^ + m/2) 7 ^(l - 7 5 )] (B.10) 

By using the following properties of the 7 -matrices: i) Tr[ 7 a 7 ^ 7 7 ..] = 0 (with [ 7 a 7 ^ 7 7 --]) an 
odd number of 7 matrices; ii) Tr[ 7 5 7 a 7 ^] =Tr[ 7 5 7 “ 7 ^ 7 7 ] = 0 ; iii) [ 7 ^, 7 5 ]+ = 0 (with [] + 
the anitconmutator); iv) ( 7 s ) 2 = 1. Then, equation (B.10) reads 

A al9 = 2p , lA Pi 7 Tr[7 A 7 a 7 7 7 /? - 7 5 7 A 7 Q 7 7 7 /3 ], (B.ll) 

where p' = 7 X p'\ x and p = 7 7 pi 7 - So, by using the following properties v) Tr[ 7 A 7 Q 7 7 7 ^] = 
4(g A “g 7/J — g A 7 g Q/3 + g A ^g“ 7 ); vi) Tr[ 7 5 7 'b a 7 7 7 ^] = —4 ie Xa1 ^, we finally have 

A “' 3 = 8(p?pf - p7p h g a/3 + P?pf) + 8 zp' lA p l 7 e A “ 7 ^ . (B.12) 

Analogously as done for the A factor, the B factor reads (using i-vi) and recalling g a sj s = j a , 
etc) 

Bap = 8>{P2 a P2p - P2P2„g a p + P2 a P2p) + ^P'^Pl^lP , (B.13) 

where p' 2 = 7 M p 2 and p, t = ^ V p 2 v - Finally, by multiplying Eqtns. (B.12) and (B.13) we have 

A al3 B a p = M[26(p\p l 2 )(pip 2 ) - :i(p\p\)(p 2 p 2 ) + 2(p' 1 p2)(y>iy4)] , (B.14) 
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where p'^Pi^ = Pip' t (i = 1,2), and we have used the important properties of the Levi-Civita 
symbols (or alternating symbols), vii) e Xai ^e\ ai p = —24; and viii) = ±1, with the + 

(—) sign is for (A, tt) and even (odd) permutation of (0,1,2,3), and vanishes if two or 
more indices are the same. This last property was used to show that the crossing-products 
in A (B.12) times B (B.13), are 0, and only two terms survive leading to (B.13). 

We proceed with the calculation of the factor-elements E and F, and notice that 
E a/3 = f A a P, and F af 3 B a p. Therefore, it leads straightforwardly to 

E 0 * = 8 (p'M - p^pi^ + pfp?) + 8 ip' 2x p h e Xa ^ (B.15) 

Fap = KPl a P2p ~ PlP2„gap + P2,Jhp) + ^PlPl^fi , 

leading directly to (after using vii-viii) 

E aP F a p = 64[26(pip' 2 )(pip 2 ) - 3(pip2)(pu4) + 2 (PiPi)(P 2 p' 2 )] • (B.16) 

With the results obtained up to here from equations (B.14) and (B.16), plus the use of 
the v s - u s elastic scattering kinematics detailed in the last section below, we finally have for 
Faa and X pp 


X aa = (—') , 2 4 1 4( nl 4 [ 26 (PlP 2) 2 - ^(PlPlf + ^PYP'if] 

\mv J 4 z m 4 cos 4 (#(y) 

X bb = ( — ') 42 4 1 4(01 J 26 (PlP 2) 2 - 3(p 1/4) 2 + 2(jpip' l ) 2 ] . 
\my J 4 z m 4 cos 4 ( 6 qy) 


(B.17) 

(B.18) 


We now calculate the last term X a b, which can be written as traces of products of up 
to nine 7 matrices (after having used all the properties i-vi) 


Xab = (C/m 4 ) 8 Trty' 1 'y a p 1 'yPp , 2 'y a p' 2 'yi 3 - ■ (B.19) 

By the use of the following properties in equation (B.19): ix) Tr[{7V] =Tr [VU] (with U,V 
4x4 matrices), and the following contraction properties of 7 -nratrices x) 7 a7°7^7 7 7 A = 
— 27 7 7^7 a ; xi) 7 a-^$7^ = 4 PQ (with P and Q four-vectors) ; it can be easily shown the 
following result (coming from the first trace term in (B.19) only, while the second trace term 
containing 7 ® vanishes by the use of the extra-property ii) 

= -(C/mVfeft W/ 2 ] = - (^) 4m 4 co 1 5 4 (aW) (P'^> 2 > 

where for the right hand side in the above equation we have used the property xii) Trpy 0 ^] = 
4 g a/3 , and the kinematical properties given in the last section. 

We note that X a f, is real, and therefore we have X* b = X a b. This last result together 
with (B.18) and (B.20) implies finally for the differential cross-section (B. 6 ), the following 
equation (we have used the elastic scattering condition p( = p l5 and the assumption of 
typical initial sterile neutrino energies (E\ « E 2 = E) 


da A 
dtt ) 


CoM 


gv V 

m v ) 


4 4 7 t 2 E 2 COS 4 (#(y) 


[60(pip 2 ) 2 - (piPif - (P 1 P 2 ) 


(B.21) 
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which, by the use of the kinematics for this process in terms of the typical particle energy 
E 2 = (p 2 + m 2 ), the three-momentum p, and the scattering angle 9 (see last section for 
details), we have finally for the differential cross-section 


da 


= (slS 


dn /CoM \ m vJ 4 4 7T 2 COS 4 (@(y) 


[29 E 2 + 60p 2 + p 4 /^ 2 (30 - cos 2 (0))] 


(B.22) 


leading directly by the use of equation (B.l) the total cross-section for our v s - v s scattering 
process in terms of p and m 

4 


_tot 
a CoM ~ 


1 


= (—) 

\my ) 4 3 7r 2 cos a (6' w 


-[29m 2 + 89p 2 + 89/3 


P 


p 2 + m 2 


(B.23) 


This is the expression we use in section 3.3 [eqn. (3.11)]. In fact, for our purposes in this 
work it suffices to consider only a single Abelian vector interaction, thus ignoring any mixing 
angle in the dark sector. This implies that we can safely set 9' w = 0. 


Majorana equation, plane wave solutions and energy projection operator 

As we have shown in section 2 of the paper, the equation of motion for our right-handed 
sterile neutrinos singlets Nr\ is given by 

+ l/2m) = 0, (B.24) 

with p being (B.24) obtained after taking the variations of the corresponding 

Majorana Lagrangian density Cn r1 w.r.t Nr i, given by equation (2.3) of the paper. Equation 
(B.24) admits plane wave solutions of the form 

Nri(x) = const. (u r (p))e~ ipx (B.25) 


where u r (p) are the constant singlet four-spinors with positive energy E and momentum p. 
Now, if we define a Majorana (positive) energy projector operator of the form 


A+(p) 


f — 1/2 m 

—m 


—p + l/2m 
m 


(B.26) 


it can be directly seen that fulfills with the property of projecting out the positive energy 
solutions from a linear combination of possible plane wave states u r ( p) and u r (p), with v r a 
negative energy plane wave state (if it would exist), i.e. 


u r (p)A + (p) = u r (p); u r (p)A + (p) = 0, (B.27) 

where u r (p) = ^r(p)7° is the adjoint, and v) r the transpose spinor. Therefore, with (B.27), 
and by the use of the completeness relation fulfilled by the constant-spinors: T, r [u ra (p)u r p (p) — 
Vra{p)v r p{p)\ = Sa/3 (see also [68], appendix A.4), it follows the following relation for our 
projector operator 13 . 

2 

A + (p) = A+^(p) = ^ u ra(p)Urp(p) , (B.28) 

r =1 

which will be largely used in our cross-section calculations, allowing us to calculate the 
X aa ,X bb ,X ab terms. 

13 Notice that this was done, in complete analogy with the Dirac case, where the projector was defined by 
Aj = (ft + m)/2m, in accordance with the Dirac equation u r (p)(ft — m) = 0, and following the same properties 
shown here, among others ([68], appendix A.4) 
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Figure 6. Kinematics for a u s - v s scattering process in the CoM system. 

v s - v s elastic scattering kinematics in the CoM system 

The kinematics of an elastic scattering process in the Center of Mass system is given by the 
following equations, accordingly with figure 2 


pi = (£',pi); Pi = {E, Pi) (B.29) 

P 2 = (E, -pi); p 2 = ( E , -pi) 

where the equalities p 1 = —p 2 and pi = —pi implies we are in the Center of Mass (CoM) 
system. Moreover, as we are under the assumption of elastic scattering conditions, the 
following equality for the three-momentum holds: p = p' , with p = |p| and p' = Ip^. Finally, 
we have used the assumption that the initial particle Energies are of the same order, i.e. 
Ei = E 2 = E. Then we must have 

P\P\ = P2P2 = E 2 -p 2 cos(0) (B.30) 

P1P2 = P2P1 = E 2 +p 2 cos (6) 

P1P2 = p'ip'2 = E 2 +p 2 

C Effective four-right-handed-Majorana-neutrino interaction coupling in 
the z/MSM 

In this appendix we would like to explain why the induced four-fermion Majorana right- 
handed neutrino terms are subleading compared to the strong self-interactions considered 
above. This stems from the very weak nature of the Yukawa couplings F a j as dictated by the 
seesaw mechanism which is assumed to be in operation here [ 6 ] so as to give a mass in the 
active neutrinos, via their mixing with the sterile (Majorana) neutrinos of Mass Mj = 1, 2, 3, 
of the form: 

v 2 

m au = F 2 I —, a = e,n,r (C.l) 


- 30 





N 


N 


+ topologically 
inequivalent 
configurations 


N 


> 

> 





N 




o 


>7 

>7 



+ topologically 
inequivalent 
configurations 


Figure 7. Generic Feynman diagrams for the effective four-fermion right-handed Majorana neutrino 
coupling, in the context of i/MSM. The topologically inequivalent graph configurations are not shown 
explicitly as we are only interested in an order of magnitude estimate. Thick (thin) continuous lines 
indicate right-handed (active) Neutrinos N (v). Dashed lines indicate the Higgs (H) held, and the 
grey blobs denote the Higgs vacuum expectation value ( H) ^ 0 . Upper pannel: One-loop, cut-off 
dependent diagram. Lower pannel: Tree-level diagram due to Majornana right-handed external line 
transmulation into an active (light) neutrino and Higgs held disappearance into the vacuum. 


where v ~ 175 GeV is the Higgs v.e.v. and the mass of the heavy neutrinos is assumed 
much larger than the induced active neutrino mass m aiy . In the z^MSM one has the mass 
hierarchy [ 2 ] 

Mi = 0(10 - 50 keV) < M 2 = M 3 = 0(1 GeV) 

and thus, assuming as a representative range (upper bound) of active neutrino masses the 
order of magnitude implied by the atmospheric experiments (which is in agreement also with 
the cosmological bound), 

ml ~ | Am atm | 2 = (2.40 + 0.12 - 0.11) x 10 -3 eV 2 , 
we obtain from (C.l) [4]: 

Fa 1 ~ 10 10 , F a 2 ; 3 ~ 1 (T 7 , (C. 2 ) 

thereby pointing towards very weak Yukawa couplings as stated previously. It should be noted 
though that it suffices to generate masses for three active neutrinos via (C.l) by assuming 
that the only non-vanishing Yukawa couplings are F a 2 3 , in which case the mixing of the 
lightest neutrino N± with the matter sector is suppressed. 

We next notice that the couplings of the Majorana neutrinos to the observable sector 
would itself yield an effective four Majorana neutrino contact interactions, via two kind of 
processes: 

• (i) through an Ultraviolet-momentum-cut-off-A dependent one particle irreducible box 
(one-loop) diagram (to leading contribution in perturbation expansion) with an active 
neutrino (assumed massless for all practical purposes) and a Higgs as internal lines, 
indicated in the upper panel of figure 7. Concentrating on the lightest right-handed 
neutrino for concreteness (thus assuming F a \ 7 ^ 0 (cf. (C.2)), the corresponding four 
- N\ -neutrino scattering amplitude that would generate the contact interaction in an 
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effective field theory framework would be proportional to 


(F a i) 4 x J 


d A k 11 1 1 

(2vr ) 4 J + 1 (p2 + k) 2 - Mjj (pi - k) 2 - M 2 U 


(C.3) 


where polarization spinors have been omitted for brevity, p\p denote incident Majorana 
neutrino external momenta, whilst p 3,4 the outgoing ones, and Mu is the mass of 
the Higgs field (assumed massive in the galactic era we are interested in, since the 
electroweak symmetry is broken at that epoch. 


• (ii) processes indicated generically in the lower panel of figure. 7, which involve the 
transmutation of an incident (or outgoing) Right-handed Majorana neutrino line into 
an active neutrino one, via a disapperance of a Higgs line into the vacuum. Such 
graphs contribute an effective four-Majorana Neutrino, whose amplitude is of order 
(up to numerical factors counting the various inequivalent graphs and ignoring spinor 
polarizations for brevity) 

V 2 

4_^_ 

[pi - P3 ) 2 - M'h 

where v is the Higgs field vacuum expectation value. 



The effective upper bound in energies of the Majorana neutrinos in the galactic centre 
and haloes cannot exceed the MeV scale (usually is at the keV. Hence, by using as an 
ultraviolet cutoff A -C Mu ~ 0(100) GeV, and taking into account that the dominant 
contributions to the amplitude (C.3) come from the region of integration pi(i = 1 ,.. .4) <C 
k <C Mu, one may estimate the total total four fermion amplitude, stemming from the sum 
of the terms (C.4), (C.3), as proportional to 


p4 1 w f d A k 11 4 A 2 

al Mfj X J (2vr) 4 ^ OC Ql M 4 


(C.5) 


which yields the effective four fermion coupling G^f for this induced interaction. With the 
value (C.2) of the Yukawa coupling F a 1 , and a cut-off A in the ball park of a few MeV, the 
effective coupling G 4 / is much smaller than the postulated strong four fermion coupling in 
our self-interacting model, to be discussed below 11 . Hence the conventional I'MSM model, 
without the inclusion of relatively strong self interactions, as proposed here, cannot provide 
an effective description of the galactic structure. 


14 Actually, in view of the smallness of F a 1 , the same conclusion holds for A even much larger than the 
electroweak scale Mw = 0(100) GeV. 
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